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Abstract 

The  evaluation  of  transmission  efficient/  of  structural  junctions  forms  an  important  part 
in  the  study  of  structure-borne  noise  as  it  provides  the  basis  for  identifying  and 
quantifying  the  vibration  paths  in  the  structure.  In  this  report,  analytical  methods  for 
evaluating  the  transmission  efficiency  of  structural  junctions  including  plate-plate  and 
plate-beam  junction  are  described.  The  calculation  of  coupling  loss  factor  from  the 
transmission  efficiency  of  a  junction  for  Statisticul  Energy  Analysis  (SEA)  is  also 
described.  Sample  calculations  of  transmission  ^dencies  on  a  number  of  structural 
junctions  are  presented.  It  is  found  that  for  typiail  tumal  ship  constructions  that  consist 
of  plates  coupled  to  light  thin  beams,  the  elastic  vibrations  of  the  beams  have  a  significant 
effect  on  the  transmission  efficiency. 
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The  Evaluation  of  Transmission 
Efficiency  and  Coupling  Loss  Factor  of 
Structural  Junctions 


1.  Introduction 

One  of  the  limiting  factors  in  the  analysis  of  vibration  transmission  through 
complex  structures,  such  as  naval  ships,  is  the  transmission  across  junctions  and 
discontinuities  in  the  structure  where  vibration  waves  are  partially  reflected  and 
partially  transmitted.  The  wave  transmission  properties  of  a  structural 
discontinuity  may  be  characterised  by  the  transmission  efficiency  which  is  defined 
as  the  ratio  between  the  transmitted  wave  power  and  the  incident  wave  power. 
The  transmission  efficiency  is  an  important  parameter  in  the  study  of  structure- 
borne  noise  since  it  provides  the  basis  to  identify  and  quantify  the  vibration  paths 
in  the  structure.  This  information  enables  appropriate  vibration  control 
techniques  to  be  applied.  The  transmission  efficiency  may  also  be  used  to 
calculate  the  coupling  loss  factor  for  Statistical  Energy  Analysis  (SEA).  SEA  is  a 
very  useful  technique  for  analysing  the  average  vibration  levels  of  complex 
interconnecting  elements,  such  as  ship  structures,  especially  at  high  frequencies. 

One  of  the  early  attempts  to  evaluate  the  transmission  efficiency  was  carried  out 
by  Cremer  (1948).  His  work  included  right-angled  plate  junctions  subjected  to 
oblique  incident  bending  waves.  Other  authors  (Wohle  et  al.,  1981;  Craven  and 
Gibl»,  1981;  and  Langley  and  Heron,  1990)  have  extended  the  analysis  to  include 
longitudinal  and  transverse  shear  waves.  The  contribution  of  these  vibration 
waves  to  structure-borne  noise  has  been  investigated  by  Lyon  (1986).  In  the 
previous  studies  involving  transmission  efficiency  of  plate-beam  structiues,  the 
stiffening  beam  was  modelled  by  using  conventional  beam  theory  and  the  effect  of 
beam  vibration  was  neglected.  While  this  so  called  l>locking  mass"  approach 
might  be  valid  for  thick  heavy  beams,  there  are  situations  where  elastic  vibrations 
of  the  beam  have  to  be  considered.  For  example,  in  plate-beam  structures  typical 
of  ship  constructions,  the  web  thickness  may  be  of  the  same  order  as  the  plate  and 
hence  the  effect  of  web  vibration  has  to  be  accounted  for  in  the  evaluation  of 
transmission  efficiency. 


In  this  report,  analytical  methods  for  evaluating  the  transmission  efficiency  of 
structural  junctions  typical  of  those  found  in  ship  constructions  are  described. 

The  effect  of  elastic  vibrations  of  the  beam  in  a  plate-beam  junction  is  investigated 
by  modelling  the  beam  as  a  finite  plate  coupled  to  a  system  of  semi-infinite 
plates.  Samples  of  calculations  for  a  number  of  junctions  are  presented. 


2.  Modelling  the  Structural  Junctions 

2.1  Mathematical  Expression  for  Vibration  Waves 

Figure  1  shows  a  schematic  diagram  of  plate-plate  junction  which  consists  of  n 
plates  coupled  along  a  line.  The  plates  are  assumed  to  be  infinite  along  the  y  and 
X,  to  x^  directions.  It  is  further  assumed  that  the  plates  are  thin  so  that  the 
boundary  conditions  can  be  applied  on  the  plate  centreline.  Plate  1  is  subjected  to 
an  oblique  incident  wave  which  can  be  either  bending  (B),  longitudinal  (L)  or 
transverse  shear  (T).  The  incident  wave  is  partially  reflected  and  partially 
transmitted  at  the  junction  as  bending,  longitudinal  and  transverse  shear  waves  as 
shown.  To  study  the  mechanism  of  wave  transmission  at  a  junction,  it  is 
convenient,  as  a  first  step,  to  consider  the  elastic  deformations  due  to  these 
reflected  or  transmitted  waves  in  an  arbitrary  plate.  Figure  2  shows  the  plate 
deformations  u,  v  and  w  along  a  set  of  local  co-ordinates  x,  y  and  z  respectively. 
Using  thin  isotropic  plate  theory,  the  following  governing  equations  of  motion  for 
plate  deformations  may  be  derived  (see  Love,  1927,  p496). 

The  bending  wave  equation  being; 

V*w  +  [l2p(l-ii^)/Eh^]d^w/dt^  =0,  (1) 

and  the  in-plane  wave  equations  are: 

d^u/dx^  +[(\-ii)/2]d^u/dy^  +[(l  +  p)/2]d^v/dxdy-lp{l-p^)/E]d^v/dt^  =0, 

(2) 

d^v/dy^  +[(l-p)/213^o/ax^  +I(l+M)/2]a^M/dx^-[p(l-H^)/£]a^M/df^  =0, 

(3) 


where  =19* /ax* 

p  =  material  density, 

4  =  Poisson's  ratio, 
h  =  plate  thickness, 

E  =  Young's  modulus. 


waves 


B  -  Bending 
L  -  Longitudinal 
T '  Transverse  shear 


Figure  2:  Co-ordinate  system  of  an  arbitrary  plate  showing  transmitted/r^ected  waves, 
plate  displacements  and  junction  forces  and  moment. 


The  in-plane  wave  equations  are  functions  of  the  plate  deformations  u  and  v.  To 
obtain  a  solution  for  these  equations,  one  can  make  use  of  the  velocity  potential  ({) 
and  stream  function  y  defined  as  follows  (a  detailed  discussion  on  the  use  of 
velocity  potential  and  stream  function  to  analyse  vibration  waves  is  given  by 
Cremer  et  al.,  1988,  p  138): 


M  =  -3^/3x-3y/3y,  (4) 

t;  =  -3(|i/3y  +  dy/3x.  (5) 

Using  equations  (4)  and  (5),  each  of  the  in-plane  wave  equations  is  reduced  to  a 
function  of  one  dependent  variable  (<|i  or  y)  only; 

v^4.-[p(i-p^)/E]a^4)/a/^=o,  (6) 

vV-I2p(i+p)/Eia^v/af^  =0.  (7) 

The  general  solutions  to  equations  (1),  (6)  and  (7)  may  be  expressed  as: 

w=\N  exp  X  +  kg^y+jait).  (8) 

<t»  =  <I>  exp  +  ;(Dt),  (9) 

y  =  ‘H  exp  (kj^  x  +  k^^y  +  ju>t).  (10) 


The  velocity  potential  ^  is  associated  with  longitudinal  waves  while  the  stream 
function  y  with  transverse  shear  waves.  For  the  solutions  to  be  valid,  the 
following  conditions  must  be  satisfied  (these  conditions  can  be  obtained  by 
substituting  equations  (8),  (9)  and  (10)  into  equations  (1),  (6)  and  (7)  respectively): 

=  ±ll2paj^  (1  )/ £//'  1^2  =  ±k^ ,  (11) 

-(kL+*Ly)  =  [p®'(l-l^')/£l  =  *^  (12) 

-(*rx+)^Ty)  =  l2p(0*(l  +  p)/El  =  k^,  (13) 

where  kg,  kj^  and  kj-  are  the  bending,  longitudinal  and  transverse  shear  wave 
numbers  respectively.  Snell's  law  states  that  the  trace  velocity  of  all  wave  types  at 
the  junction  must  be  the  same.  This  implies  that  the  y  -  component  of  wave 
numbers  (i.e.  k,^  kj^^  and  kj^)  for  the  reflected  and  transmitted  waves  of  all  plates 
must  be  the  same  as  that  of  the  incident  wave.  The  x  -  component  of  wave 
numbers  may  be  determined  from  equations  (11)  -  (13).  For  bending  waves, 
equation  (11)  yields  four  roots,  the  negative  imaginary  root  and  the  negative  real 
root  must  be  selected  since  they  represent  propagating  and  decaying  waves 
respectively  in  the  positive  x  -  direction  (i.e.  away  from  the  junction).  Similarly, 
the  solutions  for  longitudinal  and  transverse  shear  waves  must  be  negative 
imaginary.  Equations  (4),  (5)  and  (8)-(10)  give  the  elastic  deformations  of  the 
plate,  u,  V  and  w.  By  expressing  the  x  and  y  •  components  of  wave  numbers  in 
terms  of  the  reflection/transmission  angles  y,,  >1^  and  the  plate  deformations 

may  be  written  as: 
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M  =  I;\4»cosyL  exp^-jk^^xcosy^^)-ik■J.V^yJ.  exp{- jk-j-x cosy j  )1 
exp(<tyy+;o>t). 

(14) 

o  =  I-;\«&sin 7 L  exp(-  jk^x cos y l  ) -  jk^'Vcosyj.  exp(- jkjX cos y ^ )] 
exp(kyy  +  ;<Bt), 

(15) 

a;  =  { JV,  exp(  -ykgX  cos  y  g )  +  V\f^  exp(  -kg  x  7(1  +  sin*  yj)])  exp(kyy  +  ;(0t),  (16) 

where  ity  is  the  y  -  component  of  incident  wave  number.  The  last  exponential 
factor  on  the  right  hand  side  of  the  above  equations  represents  the  y  -  direction 
dependency  and  time  dependency  of  the  displacement  amplitudes.  These  factors 
are  exactly  the  same  for  all  wave  types  and  will  be  omitted  in  the  subsequent 
expressions  for  simplicity.  The  reflection/transnussion  angles  may  be  expressed 
in  terms  of  the  incident  angle  using  the  Snell's  law  as  discussed  earlier: 

kg  sinYg  =k^^  siny^  =kj  siny^.  =  ksina,  (17) 

where  a  is  the  incident  angle  and  k  is  the  incident  wave  number.  Note  that  the 
incident  wave  can  be  either  bending,  longitudinal  or  transverse  shear.  Equations 
(14) '  (16)  can  be  further  simplified  by  introducing  the  complex  wave 
displacement  amplitudes  defined  as  follows: 

4bn  =  ^2'  (18)'  (19) 

(20),  (21) 

^g  and  represents  the  complex  wave  amplitudes  of  the  travelling  and  decaying 
bending  waves  respectively  while  and  ^j.  denote  the  amplitudes  of  the 

longitudinal  and  transverse  shear  waves.  Using  equations  (17)  -  (21),  the  plate 


deformations  may  be  expressed  as: 

“  =  /k)lexp(-/kRj^x)-^r[sina/(kj.  / k)]exp(-)kJij.x),  (22) 

v  =  -^Jsma/(ki^/k)]exp(-jkRi^x)-^^lRj./(kj. /k)Jexp(-fkR^x),  (23) 

“'  =  4Bexp(-/kRgX)  +  ^gg,exp(-kRg^x),  (24) 

where  Rp  =  V((kpA)^-sin*al,  (25) 

D  =  B,  LorT  depending  on  the  wave  type, 

RgN  =  /*)*  +  sin^al  (26) 


In  equation  (25),  if  the  quantity  inside  the  square  root  is  negative,  the  exponential 
term  representing  the  travelling  wave  in  the  wave  equations  becomes  a  real 
quantity  and  wave  propagation  cannot  exist.  One  must  then  replace  the  quantity 
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>/l(jfcpA)^-sin^<*l  by  -/V[sin*o  -(kp/k)^]  in  the  solution  of  wave  amplitudes.  For  plate 
1,  the  deformations  must  also  include  the  components  of  incident  wave.  The 
elastic  deformations  due  to  the  incident  wave  may  be  obtained  in  a  similar 
manner  as  that  of  the  reflected/ transmitted  waves.  The  only  difference  being  that 
the  positive  imaginary  solution  of  the  wave  numbers  in  equations  (11)  -  (13)  must 
be  selected  since  the  incident  wave  propagates  towards  the  junction.  For  the 
purpose  of  evaluating  the  transmission  efficiency/  the  mcident  wave  may  be 
considered  as  having  a  imit  amplitude.  Hence,  for  an  incident  bending  wave,  the 


elastic  deformation  is: 

a;®  =  exp  (/ fcgX  cos  a),  (27) 

for  a  longitudinal  wave, 

u‘-  =  -  cos  a  [exp  (;\x  cosa)l,  (28) 

r/-  =  -  sin  a  (exp  (/\  x  cos  a)],  (29) 

and  for  a  transverse  shear  wave, 

uT  =  -  sin  a  [exp  (/'  k-^  x  cos  a)],  (30) 

=  cos  a  [exp  (/ itj  X  cos  a)].  (31) 


Superscripts  B,  L  and  T  denote  the  type  of  incident  wave.  If  follows  from  the 
above  analysis  that  the  elastic  deformations  of  each  plate  in  the  junction  are 
expressed  in  terms  of  four  unknowns  representing  the  complex  wave  amplitudes, 
namely,  the  amplitudes  for  longitudinal  and  transverse  shear  waves,  as  well  as 
the  travelling  and  decaying  bending  waves.  Hence,  in  a  junction  that  consists  of  n 
coupled  plates,  there  are  4n  unknowns  to  describe  the  wave  motion.  These 
unknowns  may  be  solved  by  the  appropriate  boundary  conditions. 


2.2  Boundary  conditions 
2.2.1  Plate-plate  junction 

To  consider  the  boundary  conditions  at  a  junction,  it  is  convenient  to  introduce 
the  subscript  t  to  denote  the  plate  number  (i  =  1, 2,...n).  The  compatibility  of  plate 
motions  requires  that  the  displacement  components  of  all  plates  along  a  set  of 
reference  co-ordinates  (e.g.  x,,  y  and  z,)  at  the  junction  must  be  the  same.  In 
addition,  the  rotation  about  the  y  axis  of  all  plates  should  be  equal.  The  plate 
rotation  is  given  by: 

e,  =  dw/dx^.  (32) 

The  displacement  components  of  all  plates  (u^  and  wj)  due  to  bending, 
longitudinal  and  transverse  shear  waves  are  given  by  equations  (22)  -  (24).  For 
plate  1,  the  displacements  must  also  include  components  of  the  appropriate 
incident  wave  as  given  by  equations  (27)  •  (31).  Resolving  the  displacements  of 


plate  1  along  co-ordinates  Xj,  y  and  2,  leads  to  the  following  compatibility 
equation: 


n 


“i 

cosp,.,  Q  -sinp,.  0 

“1 

0 

0 

0 

sinP,  ,  0  ,  cosP,.  ,  0 

0. 

0 

0 

0 

01 

1  =  2 

where  P,  =  angle  between  plate  i  and  plate  1. 

For  thin  isotropic  plates,  the  forces  and  moment  per  unit  length  of  plate  (along 
the  y  -  direction)  may  be  expressed  as  (see  Leissa,  1969,  p  336): 


=-[£i^  /(l-p,^)PMi /Bx.  +p,ay,. /dyl  (34) 

Fj,,  =-[£A  /2(1  +  P,)1[3m,  /By  +  dv. /dx-],  (35) 

F«  =l£A'/12(l-ji.2)][33t,;./ax.3-|.(2-p.)a^a;./ax,dy2l,  (36) 

M.  =-[£(/i/  /12(l-p.^)113^a;. /Bx.^  +n,B^w,  /By^].  (37) 


Note  that  equation  (36)  represents  the  effective  edge  force  which  consists  of  the 
shear  force  plus  the  effect  of  twisting  moment  (see  Leissa,  1969,  p  338  and  Cremer 
et  al.,  1988,  p  435  for  details).  Using  equations  (34)  -  (37),  the  forces  and  moments 
of  the  plates  can  be  obtained.  Equilibrium  requirements  for  the  balance  of  forces 
and  moments  at  the  plate  junction  lead  to  the  following  equation: 


n  n 

£xi  +  X  F«  P,.  +  ^  F„  sin  P,  =  0, 

(38) 

1=2  i=2 

P,-o 

l»] 

(39) 

Fxi  -XFx,  +X£x.  cosP,.  =  0, 

(40) 

i*2  fs2 

£m,=o. 

(41) 

<-] 


Equations  (33)  and  (38)  -  (41)  provide  the  necessary  boundary  conditions  for  the 
solution  of  wave  amplitudes.  These  equations  are  evaluated  at  the  origin  of  the 
co-ordinate  system  (i.e.,  Xj,  x^,  Xy...=  0).  A  standard  computer  routine  which 
handles  a  system  of  linear  equations  with  complex  coefficients  may  be  used. 


2.2.2  Plate-beam  junction  (thick  beam) 

Figure  3  shows  a  schematic  diagram  of  an  infinite  plate  coupled  to  a  thick  uniform 
beam.  The  thick  beam  assumption  in  the  context  of  this  report  implies  that  no 
incident  wave  is  transmitted  to  the  beam.  Also,  it  is  assumed  that  the  beam 
centroid  coincides  with  the  shear  centre  and  is  a  principal  axis.  It  is  further 

assiuned  that  the  boundary  conditions  can  be  applied  to  the  origin  of  the  co¬ 
ordinate  system  ( the  validity  of  this  assumption  is  discussed  in  section  5.2). 

Thus,  the  compatibility  requirements  for  plate  motions  in  this  case  are  exactly  the 
same  as  for  a  plate-plate  junction  (equation  (33)).  However,  the  force  and  moment 
balance  equations  ((38)  -(41))  must  be  modified  to  allow  for  the  torsional,  bending 
and  inertia  effects  of  the  beam.  For  beams  of  high  slenderness  ratio  (defined  as 
the  ratio  between  the  length  and  radius  of  gyration),  the  effects  of  shear 
deformation  and  rotary  inertia  may  be  neglected  (these  effects  have  been 
investigated  by  Langley  and  Heron,  1990)  and  the  forces  and  moment  balance 
equations  may  be  derived  by  the  conventional  beam  theory. 
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Figure  3:  A  plate  -  thick  beam  junction  showing  the  displacements  and  junction  forces 
and  moments. 


At  the  junction,  the  sum  of  the  forces  and  moments  due  to  all  plates  may  be 
resolved  along  the  beam  co-ordinates  Xy  y,  and  expressed  as  the  beam  forces 
and  moment  (see  Figure  3): 

-f»i  sinp>  cos(P,.  -P»)+^f„  sin(P,  -p^),  (42) 

i-2  i>2 


(43) 


II  m 

-F,,  sinPj  +  f,,co*Pj-5^f,,  sin(P, cos(P, -p^),  (44) 

i«>2  f«2 

(45) 

The  motions  about  the  beam  centroid  and  6^  may  be  expressed  in  terms 

of  the  elastic  deformations  of  plate  1; 


s  M,  COS  p>-  w,  sin  P>, 

(46) 

V^  =  Vj, 

(47) 

Wy  =  u,  sin  Pj  +  u>,  cos  P^  +  6,S, 

(48) 

(49) 

where  Pj  =  angle  between  the  beam  and  plate  1, 

S  =  distance  between  beam  centroid  and  the  junction. 

The  equilibrium  of  forces  and  moments  at  the  junction  must  allow  for  the 
torsional,  bending  and  inertia  effects  of  the  beam.  Consider  the  balance  of  forces 
in  the  Xj  direction,  the  beam  force  is  augmented  by  the  shear  force  as  a  result  of 
beam  bending  in  the  Xj-y  plane  (a  similar  argument  exists  for  forces  in  the  Zj, 
direction).  Summation  of  forces  along  the  beam  co-ordinates  leads  to  the 
following  force  balance  equations: 

(50) 

(51) 

(52) 

where  = 

mass  per  unit  length  of  beam. 

£»  = 

Young's  modulus  of  beam. 

/)*  = 

second  moment  of  area  of  beam  about  axis  x^ 

= 

second  moment  of  area  of  beam  about  an  axis  parallel  to  and 
passes  through  the  beam  centroid. 

The  variation  in  plate  rotation  8^  along  the  y  -  axis  causes  the  beam  to  twist  and 
results  in  a  torsional  moment.  Consider  the  equilibrium  of  moments  about  a  line 
parallel  to  the  y  •  axis  and  passes  through  the  ^am  centroid: 

-M,*F^S+T,d\/dy^  (53) 
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where  =  torsional  stiffness  of  the  beam. 

Ii,  =  moment  of  inertia  per  unit  length  about  beam  centroid. 

Equations  (33)  and  (50)  -  (53)  represent  the  boundary  conditions  for  a  system  of 
plates  coupled  with  a  thick  beam.  Again,  these  are  sufficient  to  permit  solution 
for  the  wave  amplitudes  of  interest. 


2.3  Mathematical  Model  for  a  Thin  Beam 

Some  engineering  structures  (e.g.  ships,  aircraft)  quite  often  involve  the  use  of 
thin  beams  to  reinforce  plate  elements.  A  thin  beam  in  the  context  of  this  paper 
implies  that  the  beam  thickness  is  of  the  same  order  as  that  of  the  plate  element 
and  is  therefore  subjected  to  bending  and  in-plane  waves.  A  schematic  diagram 
of  the  structure  is  shown  in  Figure  4.  The  analysis  of  this  type  of  structures  may 
be  carried  out  by  assuming  that  the  thin  beam  behaves  as  a  Hnite  plate  with  waves 
travelling  in  both  the  positive  and  negative  x  -  directions.  Referring  back  to 
equations  (11)  -  (13),  the  solutions  to  wave  motions  of  the  finite  plate  in  this  case 
must  include  the  positive  and  negative  roots.  This  results  in  eight  unknown 
complex  wave  amplitudes  (instead  of  four  unknowns  as  in  the  case  of  an  infinite 
plate).  By  denoting  the  web  as  the  second  plate  in  a  plate-plate  junction,  the 
expressions  for  bending  and  in-plane  waves  may  be  derived  in  a  procedure 
similar  to  those  described  in  section  2.1: 


“2  =  I  /  ^ ) il ^i.2  exp( .^2 ) -  ^  £.2  exp ( j jr^ ) ] - 

(sina/fkj-j  /  expf-jkR^x^  )  +  ^ >2  expf/lcR^jXj )], 

i/j  =  -(sina/(fcL2 exp(-7M^2^2  t-z  expf/M^jXj  ))- 
l^r2  /(*^r2  /^)11^T2  eXp( ) J, 


(54) 


(55) 


?<’,  exp(-yMg,.v,)  +  ^pj^,3  P’‘P(-*^e.\2-''2  )  +  '’’‘P(y*Vv''2  + 

4  8.V2  exp(kRg.^.,.v,). 


(56) 


The  additional  four  unknowns  in  a  plate  -  thin  beam  junction  %'b2,  ^bnz’ 

^^2  represent  waves  that  travel  in  the  negative  Xj  -  direction  of  the  finite  plate. 
Four  additional  boundary  conditions  are  thus  required  to  solve  the  wave  motions. 
These  boundary  conditions  may  be  obtained  by  considering  the  force  and  moment 
balance  at  the  end  of  the  finite  plate.  Substituting  equations  (54)  -  (56)  into  (34)  - 
(37)  and  evaluating  these  equations  at  x,=  L,„  gives  the  forces  and  moment.  For 
the  structural  junction  shown  in  Figure  4,  the  forces  and  moment  must  vanish  at 
the  free  end  of  the  plate. 

Figure  5  shows  a  typical  reinforced  plate  structure  used  in  naval  ship 
constructions.  The  stiffening  beam  in  this  case  consists  of  a  web  and  flange.  For 
the  analysis  of  vibration  transmission,  the  flange  may  be  considered  as  another 
finite  plate  attached  to  the  web  and  the  flange  with  the  web  junction  analysed 
using  a  similar  procedure  as  previously  described.  Alternatively,  if  the  flange  is 
thick  compared  with  the  web  (say,  or  the  order  of  twice  the  web  thickness),  one 
may  assume  that  the  flange  behaves  as  a  thick  beam  attached  to  the  web  and  the 
analysis  carried  out  in  a  procedure  similar  to  those  described  in  section  2.2.2. 
Figure  6  shows  the  forces  and  moments  at  the  flange  -  web  junction.  The  force 
and  moment  balance  equations  (evaluated  at  Xj  =  L^)  may  be  expressed  as: 
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F^j=mfd\/dt\  (57) 

F^2=mfd\/dt\  (58) 

F,j  -dF^/dy  =  / ^/^  (59) 

Mj+dMf/dy  =  Jfd%  /  3/^ ,  (60) 

where  subscript  /represents  the  flange,  and 

dF^/dy  =  Efl^d*W2  /dy*, 
dM^/dy  =  T^d%/dy^. 


Equations  (57)  -  (60)  represents  the  additional  boundary  conditions  for  the 
solution  of  wave  amplitudes. 

Z. 


1 


Figure  4;  Schematic  diagram  of  a  plate  -  thin  beam  junction. 


Figures:  Periodically  stiffened  panel. 
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3.  Transmission  and  Reflection  Efficiencies 

The  above  analysis  solves  the  wave  amplitudes  of  plate  •  plate  and  plate  •  beam 
junctions  and  leads  to  the  calculation  of  wave  power.  The  wave  power  per  unit 
length  of  a  junction  subjected  to  normal  waves  (perpendicular  to  the  junction) 
may  be  expressed  as  the  energy  per  unit  area  multiplied  by  the  group  velocity  (see 
Cremer  et  al,  1988  p  109).  For  oblique  waves,  the  expression  must  be  modified  by 
multiplying  it  with  the  cosine  of  the  wave  angle  since  it  is  the  projected  length  of 
the  boundary  line  that  effectively  intercepts  the  wave.  Figure  7  shows  the 
interception  of  an  oblique  incident  wave  for  a  plate  -  plate  junction.  Recall  that 
the  incident  wave  has  a  unit  amplitude,  the  power  due  to  an  incident  wave  at  an 
angle  a  is  given  by: 


=  m,  o)^c^,  cos  o.  (61) 

Similarly,  the  transmitted/reflected  power  may  be  expressed  as: 

Ppi  =  I  5d,  I  cos7d,.  (62) 

The  transmission/reflection  angle  Yq,  may  be  expressed  in  terms  of  the  incident 
angle  a  using  Snell’s  law  (equation  (17)).  The  power  expression  thus  becomes; 

^Di  (63) 

where  =  mass  per  unit  area  of  plate  i, 

D  =  B,LorT  depending  on  the  wave  type, 

=  group  velocity  of  plate  i, 

=  2  Cg^  for  bending  waves, 

=  Cu  for  longitudinal  waves, 

=  Cji  for  transverse  shear  waves. 
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Incident  wave 


Junction 


Figure  7:  Interception  of  an  oblique  incident  wave  at  a  junction. 

Note  that  for  bending  waves,  the  transmitted  or  reflected  power  in  a  plate  is 
determined  by  the  far  field  travelling  waves  only.  The  near  field  waves  decay 
exponentially  along  the  direction  x,  and  carry  no  time  averaged  power.  Also,  no 
power  is  transmitted  by  a  finite  plate. 

The  transmission/reflection  efficiency  is  defined  as  the  ratio  of  the 
transmitted/reflected  wave  power  to  the  incident  wave  power  and  is  a  function  of 
the  incident  angle  a: 

rf'  (a)  =  Ppi  / 

where  q  and  r  represent  the  wave  type  of  the  incident  and  generated  waves 
respectively,  and  i  represents  the  carrier  plate  of  the  generated  waves. 
Conservation  of  energy  requires  that  the  sum  of  all  transmission  and  reflection 
efficiencies  to  be  equal  to  one. 


4.  Coupling  Loss  Factor 

The  coupling  loss  factor  used  in  SEA  defines  the  amoimt  of  energy  flow  from  one 
element  to  the  other.  It  can  be  shown  (Lyon,  1975,  p  91)  that  for  two  coupled  sub¬ 
systems,  the  power  lost  by  sub-s3rstem  1  due  to  coupliitg  to  sub-system  2  is 
proportiorwl  to  the  energy  of  su^ystem  1  aiKl  may  be  expressed  in  terms  of  flte 
coupling  loss  factor  as  follows; 

P|2  (65) 

where<£,>,  *  time  averaged  erte^  of  sub-system  1, 

i),2  *  coupling  1m  factor  between  sub-systems  1  and  2. 
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Consider  an  SEA  system  that  consists  of  two  coupled  plates  as  shown  in  Figure  7. 
If  plate  1  carries  a  ^tfuse  vibratimt  field  incident  on  the  junction,  the  total  power 
transmitted  to  plate  2  can  be  obtained  by  multiplying  the  power  of  plate  1  with 
the  transmissicm  efficiency  and  then  average  the  results  over  the  entire  range  of 
incident  angles  (see  Cremer  et  al.,  1988,  p.  426): 

r«/*  _ 

^12  >/A,2k)J  ^^^T^(o)co80uio. 

=  (c„L,  <£,  >/A,*)f'Tf  {o)<f(sino), 

'Jo  *  (66) 

»(c^,L,  <E,  > /A, 

where  L,  =  coupling  line  length, 

=  mean  transmission  efficiency, 

A,  =  area  of  plate  1. 

From  equatioru  (65)  and  (66),  the  coupling  loss  factor  between  two  coupled  plates 
is  given  by: 

’ll:  (67) 

From  the  above  considerations  on  coupling  loss  factor,  it  is  evident  that  the  mean 
transmission  efficiency  is  a  useful  parameter  for  characterising  the  wave 
transmission  properties  of  a  junction  subjected  to  an  incident  diffuse  vibration 
field.  The  mean  transmission  efficiency  may  be  (rf>tained  by  a  numerical 
integration  procedure  and  incorporated  into  a  computer  program.  The  sample 
calculations  presented  in  the  following  section  are  based  on  mean  transmission 
efficiencies. 


5.  Applications  to  Ship  Structural  Junctions 


5.1  Plate-Plate  functions 

Figure  8  shows  a  schematic  diagram  of  a  section  of  ship's  structure  which  mainly 
consists  of  plate  joints.  The  mean  transmission  ef  Hciencies  for  junctions  A,  B  and 
C  were  calculated  by  the  nwthematical  model  described  in  Section  2.2.1  and  are 
plotted  in  Figures  9, 10  and  11  respectively.  In  all  tfuee  cases  studied,  it  can  be 
seen  that  the  significance  of  hvplane  vibraticm  in  wave  transmission  increases 
fvith  frequency.  For  example,  at  8000  Hz,  about  10-15%  of  the  incident  wave 
power  is  tnra^tted  by  in*plane  vibration.  Alfiwugh  in-plane  motions  are  not 
coupled  efficiently  to  the  smind  field,  ttiey  may  propagate  fiuough  the  structure 
and  transform  into  bending  ntotions  at  a  structural  discontinuity.  Hence  in-plane 
motions  nnay  be  considered  as  a  flanking  path  to  bending  motions  and,  if  that 
flanktog  pa^  is  ignored  in  the  analysis,  the  results  may  lead  to  an 
uiKlerestinuition  of  the  transmitted  structure-borne  itc^. 


Mean  transmission  efficiency 


Figure  8:  Schematic  diagram  of  a  ship's  section. 


Figures:  Wave  power  transmissim  for  junction  A. 
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Mean  transmission  etficiency  ^  Mean  transmission  efficiency 


Frequency,  Hz 


10:  Wave  power  transmission  for  junction  B. 


Frequency.  Hz 


Figure  11:  Wave  power  transmission  for  junction  C. 
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5.2  Plate-Beam  Junctions 


The  mathematical  models  for  a  thick  beam  (Section  2.2.2)  and  a  thin  beam 
(Section  2.3)  coupled  to  thin  plates  were  applied  to  the  plate-beam  junction  shown 
in  Figure  12.  Since  the  beam  thickness  chosen  in  this  example  is  the  same  as  the 
plate,  the  beam  would  vibrate  due  to  the  incident  wave  and  it  is  reasonable  to 
argue  that  the  thin  beam  model  would  give  a  more  acouate  prediction  of  the 
transmission  efficiency.  Figure  13  shows  the  bending  wave  transmission 
efficiency  of  the  junction  calculated  by  both  models.  The  thick  beam  model 
predicts  a  low-pass  characteristic  of  the  plate-beam  junction  and  underestimates 
the  transmission  efficiency  at  frequencies  above  500  Hz.  The  effect  of  resonant 
bending  frequency  of  the  thin  beam  on  wave  transmission  can  be  observed. 

As  a  second  example  on  plate-beam  junctions,  the  beam  thickness  in  Figure  10  is 
increased  to  20  mm.  Figure  14  shows  the  calculated  transmission  efficiency. 

Below  1  kHz,  the  agreement  between  the  thick  beam  model  and  the  thin  beam 
model  is  reasonable.  At  higher  frequencies,  the  thin  beam  model  predicts  a  higher 
transmission  efficiency,  possibly  due  to  the  effect  of  plate  resonance.  It  should  be 
noted  that  the  mathematical  models  used  in  this  report  are  based  on  a  thin  plate 
theory  and  the  assumption  that  the  boundary  conditions  can  be  applied  on  the 
beam/ plate  centreline.  These  assumptions  may  not  be  justified  at  high 
frequencies  where  the  cross  sectional  dimensions  of  the  junction  is  not  negligible 
compared  with  the  bending  wavelength.  Cremer  et  al.  (1988,  p  115)  suggested 
that  the  thin  plate  theory  may  be  used  if  the  bending  wavelength  is  longer  than 
six  times  the  plate  thickness.  For  the  present  example,  this  converts  to  a 

(frequency  of  approximately  26  kHz.  "nte  effect  of  plate  offset  from  the  centreline 

of  a  thick  beam  may  be  considered  by  modifying  the  compatibility  and 
i  equilibrium  equations.  This  approach  has  been  carried  out  and  reported  by 

j  Langley  and  Heron  (1990).  Dropite  the  assumptions  used  in  the  models,  the 

I  present  analysis  shows  that  the  conventional  heavy  beam  theory  may  lead  to  a 

i  serious  underestimation  of  the  transmission  efficiency  when  applied  to  thin  beam 

junctions. 


6.  Conclusions 

Analytical  methods  for  evaluating  the  transmission  efficiency  of  structural 
junctions  have  been  presented.  The  mean  transmission  efficiency  may  be  used  to 
identify  and  quantify  vibration  transmission  paths  in  a  junction.  A  study  of  the 
mean  transmission  efficiency  of  plate  junctions  shows  tluit  the  effect  of  in-plane 
motions  is  significant  in  structure-borne  noise  transmission,  especially  at  high 
frequencies.  For  structural  junctions  that  consist  of  a  thin  beam,  the  elastic 
vibrations  of  beam  plays  an  important  part  in  wave  transmission  and  should  be 
considered  in  the  analysis  of  structure-borne  noise. 
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Figure  12:  A  simple  plate  -  thin  beam  structure. 


100  1000  1 
Frequency.  Hz 

Figure  13:  Bending  wave  transmission  for  a  plate  -  thin  beam  junction. 


Frequency,  Hz 


Figure  14:  Bending  wave  transmission  for  a  plate  -  thick  beam  junction. 
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Appendix 
List  of  Symbols 


A 

B 


D 


<E,> 


‘xb 


‘zb 


i 

k 


area  of  plate  i 
bending  wave 

bending,  longitudinal  and  transverse  shear  wave 
velocities  of  plate  / 

group  velocity  of  plate  i 

subscript  to  denote  bending  (B),  longitudinal  (L)  and 
transverse  shear  wave  (T) 

elastic  modulii  of  plate  i,  beam  and  flange 

time  averaged  energy  of  sub-system  i 

component  of  internal  plate  forces  per  unit  length  in 
the  Xj,  y  and  2,  directions  of  plate  i 

thickness  of  plate  i 

second  moment  of  area  of  beam  about  axis 

second  moment  of  area  of  beam  about  an  axis  parallel 
to  2^  and  passes  through  the  beam  centroid. 

second  moment  of  area  of  flange  about  axis  Xj 

moment  of  inertia  per  unit  length  about  longitudinal 
centroidal  axis  of  beam  and  flange 

subscript  to  indicate  plate/sub-system  number 
(i  =  l,2,3 . n) 

complex  operator 

wave  number  of  incident  wave 


y-  component  of  incident  wave  number 

kg^kg^kfj  X  and  y  -  components  of  wave  numbers  for  bending, 

longitudinal  and  transverse  shear  waves 

^Ly'^x'^Ty 

bending,  longitudinal  and  transverse  shear  wave 
numbers  of  plate  f 
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L 

4 

K 

Mi 

M^Mj 

nti 

n 

^INC 

^Dl 

Pn 


longitudinal  wave 

length  of  coupling  line  at  a  junction 

width  of  web 

internal  bending  moment  per  unit  length  of  plate  i, 
the  moment  vector  is  in  the  y  -  direction 

moments  per  unit  length  about  longitudinal 
centroidal  axis  of  beam  and  flange 

mass  per  unit  length  of  beam  and  flange 

mass  per  unit  area  of  plate  t 

number  of  plates  in  a  junction 

power  of  incident  wave 

transmitted  or  reflected  wave  power  of  plate  i 

power  lost  by  sub-system  1  due  to  coupling  to  sub¬ 
system  2 


q,r 


superscripts  to  indicate  wave  type  of  incident  and 
generated  waves  respectively 


Rsi'^U’^Ti  parameters  to  .represent  the  cosine  function  of  wave 

angle  for  bending,  longitudinal  and  transverse  shear 
waves  of  plate  i  as  defined  by  equation  (25) 

Rgfjj  parameter  to  represent  the  cosine  function  of  wave 

angle  for  the  near  field  bending  wave  of  plate  i  as 
defined  by  equation  (26) 

S  distance  between  the  beam  centroid  and  the  junction 

T  transverse  shear  wave 

Tg,T^  torsional  stiffness  of  beam  and  flange 

f  time 


u,  displacement  of  plate  i  in  the  -  direction  due  to 

transmitted  or  reflected  waves 

uK  displacement  of  plate  1  in  the  x  -  direction  due  to 

incident  longitudinal  and  transverse  shear  waves 
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Vi 

v^.v^ 

a 

P. 

4* 


displacement  of  plate  i  in  the  y  -  direction  due  to 
transmitted  or  reflected  waves 

displacement  of  plate  1  in  the  y  -  direction  due  to 
incident  longitudinal  and  transverse  shear  waves 

displacement  of  plate  t  in  the  Zj  -  direction  due  to 
transmitted  or  reflected  waves 

displacement  of  plate  1  in  the  z  -  direction  due  to 
incident  bending  wave 

system  of  co-ordinates  of  plate  t 

system  of  co-ordinates  of  beam 

incident  wave  angle 

angle  between  plate  1  and  plate  > 

angle  between  plate  1  and  beam 

velocity  potential 


wave  angles  of  transmitted  /  reflected  bending, 
longitudinal  and  transverse  shear  waves 

T|i2  coupling  loss  factor  between  subsystem  1  and  2 

p,  Poisson's  ratio  of  plate  i 

0,  angular  displacement  of  plate  i  about  the  y  -  axis 

T(a)  transmission/reflection  efficiency  as  a  function  of  the 

incident  wave  angle,  the  efficiency  is  defined  as  the 
ration  of  the  transmitted/reflected  wave  power  to  the 
incident  wave  power 

mean  transmission/reflection  efficiency 

0)  circular  frequency 


^Bi  '^BNi  '^Ti  complex  Wave  displacement  amplitudes  for  travelling 

and  decaying  bending  waves,  as  well  as  longitudinal 
and  transverse  shear  waves  of  plate  i  travelling  in  the 
positive  X,  direction 

^'bi  >4' Li  same  definition  as  above  but  with  waves  travelling  in 

the  negative  -  direction 


V 


stream  function 
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